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I.  INTRODUCTION 


Application  of  the  Parabolized  Navier-Stokes  (PNS)  technique  to  predict 
the  three-dimensional  viscous  flow  about  slender  bodies  of  revolution  at 
supersonic  speeds  and  snail  angles  of _attack  has  been  performed  in  the  past 
with  a  considerable  degree  of  success.1”3  Of  particular  significance  to  the  US 
Army  is  the  capability  that  this  technique  provides  in  the  prediction  of  the 
Magnus  force  and  moment  for  spinning  shell.1’2  In  light  of  this  success,  the 
PNS  technique  has  been  extended  to  examine  the  aerodynamics  of  more  complexly 
configured  projectiles,  such  as  finned  bodies.  This  paper  documents  the 
application  of  the  PNS  technique  to  a  finned  projectile  configuration,  resemb¬ 
ling  a  long  L/D  kinetic  energy  (KE)  penetrator  projectile,  and  to  a  standard 
spinning  projectile. 

The  PNS  technique  is  a  space-marching  procedure  in  which  the  solution  is 
marched  spacially  down  the  body  in  the  main  direction  of  the  flow.  It  can  be 
applied  to  problems  for  which  the  flow  is  supersonic  and  does  not  contain 
imbedded  regions  of  subsonic  flow  or  regions  where  the  flow  separates  in  the 
marching  direction.  For  this  class  of  problems,  the  PNS  technique  can  be 
computationally  more  efficient,  requiring  less  computer  time  and  storaqe, 
compared  with  three-dimensional  Navier-Stokes  time-dependent  procedures. 

As  the  technique  has  evolved  from  the  shock  capturing  code  applied  by 
Sturek,  et  al.,1’3  significant  improvements  have  been  made,  enablinq  this 
technique  to  be  applied  to  the  finned  projectile  configuration  of  interest. 
These  improvements,  several  of  which  have  been  made  by  Rai,  Chaussee,  et 
al.,4’5  include  implementation  of  implicit  shock-fitting  boundary  conditions, 
implicit  smoothing  terms,  a  force  package  cast  in  generalized  geometry  form, 
and  an  elliptic  grid  generator  which  produces  the  qrid  over  the  finned  portion 
of  the  projectile.  The  implementation  of  the  implicit  shock-fitting  boundary 


1.  U.  B.  Sturek,  D.  C.  Mylin,  and  C.  C.  Bush,  "Computational  Parametria  Study 
of  the  Aerodynamics  of  Spinning  Bodies  at  Supersonic  Speeds,"  U.S.  Army 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Growid,  Maryland,  ARBRL- 
TR-02358,  August  1981.  (AD  A106074) 

2.  W.  B.  Sturek,  and  L.  B.  Schiff,  "Computations  of  the  Magnus  Effect  for 
Slender  Bodies  in  Supersonic  Flotf,"  U.S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02384 ,  December 
1981.  (AD  A110016 ) 

3.  L.  B.  Schiff,  and  W.  B.  Sturek,  "Numerical  Simulation  of  Steady  Supersonic 
Flow  Over  an  Ogive  Cylinder  Boattail  Body,"  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  vrn:ini  Ground,  Maryland,  ARBRL-TR-02363 ,  September 
1981.  (AD  A106060) 

4.  M.  M.  Rai,  and  D.  S.  Chaussee,  "New  Implicit  Boundary  Procedures:  Theory 
and  Applications,"  AIAA  Paper  No.  83-0123 ,  21st  Aerospace  Sciences 
Meeting,  January  1986. 

6.  M.  M.  Rai,  D.  S.  Chaussee,  Y.  M.  Risk,  "Calculation  of  Viscous  Supersonic 
Fionas  over  Finned  Bodies,"  AIAA  Paper  No.  83-1667 ,  Danvers,  MA ,  July  1983. 
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conditions  has  contributed  to  faster  rates  of  convergence  for  the  conical 
starting  procedure  and  allowed  larger  stepsizes  to  be  taken  in  the  marching 
direction  without  loss  of  accuracy  or  stability.  Use  of  implicit  smoothing 
has  been  required  in  order  to  obtain  solutions  over  the  finned  portion  of  the 
body. 

In  light  of  these  improvements,  application  of  the  technique  is  first  made 
to  a  standard  shell  configuration  (S0CBT)  to  establish  a  benchmark  for  the 
code,  which  is  henceforth  denoted  as  FINPNS.  Comparison  with  previous  PNS 
calculations  and  wind  tunnel  data  is  made  and  the  new  capabilities  evaluated. 
Results  for  the  finned  projectile  configuration  are  then  shown  and  compared 
with  the  NSWC  Euler  code  (SWINT)  and  experimental  data  base  values. 

II.  COMPUTATIONAL  TECHNIQUE 

Calculation  of  the  flow  field  over  the  body  is  accomplished  using  the 
parabolized  Navier-Stokes  technique.  The  PNS  technique  allows  the  solution  to 
be  spacially  marched  along  the  body  in  the  main  flow  direction  due  to  the 
parabolic  nature  of  the  governing  equations.  An  initial  plane  of  data  is 
required  to  begin  the  space  marching  procedure  and  may  be  obtained  either  from 
an  auxiliary  calculation  or  from  a  conical  starting  procedure,  as  has  been 
done  for  the  results  presented  here.  Both  the  space  marching  and  conical 
starting  procedures  are  outlined  below. 

A.  Space  Marching  Procedure 

The  thin-layer  Parabolized  Navier-Stokes  computational  technique  developed 
by  Schiff  and  Steger6  has  been  employed  to  calculate  the  flow  downstream  of 
the  nose.  The  governing  steady  thin-layer  equations  in  strong  conservative 
form  and  generalized  coordinates  are  written  below. 

A 

;\r  A  A  A 

s  3F  3G  l  as  m 

35  3n  ac  '  Re  H  '  1 

where  5,  n,  c  are  the  generalized  coordinate  variables  as  displayed  in  Figure 

1. 

5  =  5(x)  is  the  longitudinal  (marching)  coordinate 
n  =  n(x,.y,z)  is  the  ci rcumferential  coordinate 
C  =  c(x,.y,z)  is  the  near  normal  coordinate 

This  vector  equation  represents  the  thin-layer  approximation  to  the  equa¬ 
tions  of  mass,  momentum,  and  energy  conservation  in  the  three  coordinate 

A  A  A 

directions.  The  invisid  flux  vectors  E$,  F,  and  G  and  the  matrix  of  viscous 


L.  B.  Schiff ,  and  J.  L.  Steger ,  "Numerical  Simulation  of  Steady  Supersonic 
Viscous  Flouf'  AIAA  Paper  No.  79-0130 ,  17th  Aerospace  Sciences  Meeting , 
January  1979. 


terms,  S,  are  functions  of  the  dependent  variables  represented  by  the  vector, 
q(p,  pu,  pv,  pw,  e),  where  p  is  the  density,  u,  v,  and  w  are  the  velocity  com¬ 
ponents  in  the  three  spacial  directions  x,  y,  and  z,  and  e  is  the  total  energy 
per  unit  volume. 

The  parabolized  Navier-Stokes  equations  are  solved  us’ng  a  conservative, 
approximately  factored,  implicit,  finite-difference  numerical  algorithm  as 
formulated  by  Beam  and  Warming. 7  Further  details  of  the  numerical  method  may 
be  found  in  Reference  6.  Fitting  of  the  outer  bow  shock  has  been  performed  in 
these  calculations,  and  details  of  the  implicit  boundary  procedure  as  imple¬ 
mented  by  Rai  and  Chaussee  may  be  found  in  Reference  4.  A  fully  turbulent 
boundary  layer  has  been  simulated  in  each  of  the  reported  calculations  using  a 
two- layer  eddy  viscosity  model.8’9 


B.  Conical  Starting  Solutions 

The  initial  plane  of  data  required  to  begin  the  marching  procedure  is 
obtained  using  the  marching  code  by  assuming  conical  flow  at  the  tip  of  the 
projectile.  By  selecting  a  conical  grid  and  initially  setting  the  flow  field 
variables  to  the  free  stream  values,  the  solution  is  marched  one  step  down  the 
body.  The  solution  is  then  scaled  back  to  the  original  station  according  to 
the  conical  flow  assumption  and  again  inarched  a  single  step.  This  procedure 
is  repeated  until  a  converged  solution  is  obtained.  The  convergence  criterion 
for  the  conical  starting  solutions  applied  here  was  that  the  change  in  density 
between  successive  iterations  was  less  than  It)'5  times  the  free  stream  value 
for  each  of  the  points  on  the  body.  This  converged  solution  is  then  used  as 
the  initial  plane  of  data  in  the  marching  procedure. 

It  should  be  noted  that  for  calculations  involving  spinning  projectiles, 
the  conical  starting  procedure  introduces  a  small  error  :ince  the  circumferen¬ 
tial  velocity  at  the  body  surface  changes  with  longitudinal  position;  viola¬ 
ting  the  conical  flow  assumption.  This  error  is  small,  however,  and  the 
correct  circumferential  velocity  at  the  body  surface  is  accounted  for  as  the 
solution  is  marched  downstream. 

III.  RESULTS 

A.  SOCBT  Configuration 

In  order  to  evaluate  the  new  capabilities  of  the  code,  results  were 
obtained  for  a  secant-ogi ve-cyl inder-boattai 1  (SOCBT)  projectile  configuration 


7.  R,  Beam ,  and  R.  F.  Warming ,  "An  Implicit  Factored  Scheme  for  the 
Compressible  Navier-Strkes  Equations »"  AIAA  Journal,  V ol .  16 ,  No.  4,  1978 , 
pp.  85-129. 

8.  B.  S.  Baldwin,  and  H.  Lomax,  " Thin  Layer  Approximation  and  Algebraic  Model 
for  Separated  Turbulent  Floust"  AIAA  Paper  No.  78-257 ,  16th  Aerospace 
Sciences  Meeting,  . January  1978. 

9.  D.  Degani,  and  L.  B.  Schiff,  "Computation  of  Supersonic  Viscous  Flews 
Around  Pointed  Bodies  at  Large  Incidence >"  AIAA  Paper  No.  83-0034,  21 st 
Aerospace  Sciences  Meeting,  January  1983. 


for  which  a  significant  amount  of  wind  tunnel  and  computational  results  exist. 
A  schematic  of  this  projectile  configuration  is  shown  in  Figure  2. 

The  grid  for  the  current  FINPNS  computations  consisted  of  45  exponentially 
stretched  points  in  the  radial  direction  from  the  body  to  the  shock  and  an 
equal  spacing  of  points  circumferentially  around  the  body  at  10  degree 
intervals.  The  benchmark  PNS  calculations2  for  the  SOCBT  configuration  utili¬ 
zed  the  same  circumferential  spacing  but  used  50  points  in  the  radial  direc¬ 
tion,  clustering  half  the  points  near  the  body  with  the  other  half  equispaced 
to  a  distance  slightly  beyond  the  shock  to  enable  an  outer  shock  capturing 
technique  to  be  applied. 

Conical  starting  solutions  for  the  SOCBT  were  generated  at  a  distance  4% 
of  the  total  body  length  from  the  nosetip.  Stepsizes  for  the  conical  starting 
procedure  were  adjusted  to  be  as  large  as  possible  while  still  yielding  stable 
solutions  and  accurate  results  in  the  marching  mode.  As  noted  in  Reference  4, 
one  criterion  in  selecting  the  proper  stepsize  is  to  maintain  the  maximum 
Courant  number  in  the  flow  field  at  or  below  12.  The  maximum  Courant  number 
for  the  current  calculations  was  not  permitted  to  exceed  10.  In  some  cases, 
stepsizes  for  the  current  calculations  were  2  to  5  times  larger  than  those 
used  for  the  previous  PNS  solutions  to  which  comparisons  were  made.  (A 
primary  concern  for  the  previous  calculations  was  avoidance  of  departure  solu¬ 
tions  due  to  too  large  stepsizes.)  A  significant  decrease  in  the  number  of 
steps  for  convergence  was  observed  compared  with  the  previous  calculations, 
resulting  in  a  substantial  reduction  in  the  CPU  time  for  the  procedure  and  is 
attributable  to  the  implementation  of  the  implicit  boundary  conditions  and 
smoothing  terms. 

Results  have  been  obtained  for  Mach  numbers  2,  3,  and  4,  at  2  degree  angle 
of  attack,  with  and  without  spin  for  flow  conditions  duplicating  those  of  the 
experiments.10’11  All  SOCBT  results  were  generated  on  a  CDC  7600  computer 
with  a  speed  of  2.1  CPU  sec/step  for  the  stepback  procedure  and  4.2  CPU 
sec/step  for  the  marchinq  procedure. 

Longitudinal  surface  pressure  distributions  for  the  Mach  3  non-spin  case 
along  the  windward  and  leeward  sides  are  shown  in  Figure  3.  The  comparison  is 
good  between  the  FINPNS  computation  and  experiment,  most  notably  in  the  vicin¬ 
ity  of  the  discontinuities  in  streamwise  surface  curvature.  Figure  4  shows 
circumferential  pressure  distributions  at  four  axial  locations  on  the  projec¬ 
tile  compared  to  experiment.  The  trends  at  all  four  axial  locations  agree 
well,  and  at  the  fourth  station  (on  the  boattail)  the  maximum  disagreement 


10.  R.  P.  Reklis ,  and  W.  B.  Sturek,  " Surface  Pressure  Measurements  on  Slender 
Bodies  at  Angle  of  Attack  at  Supersonic  Speeds,"  U.S.  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground ,  Maryland ,  ARBRL-MR-02876 , 
November  1978.  (AD  A064097) 
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of  Static  and  Magnus  Aerodynamic  Coefficients  for  Projectile  Shapes  with 
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between  computation  and  experiment  is  less  than  2%.  For  this  case,  the  march¬ 
ing  stepsize  was  increased  by  5%  every  10th  step,  so  that  the  stepsize  varies 
from  .01  cal.  at  the  starting  plane  to  .037  cal.  at  the  end  of  the  projectile. 


The  effect  of  marching  stepsize  on  the  static  forces  (sign  conventions 
shown  in  Figure  5)  for  the  Mach  3  spin  case  is  shown  in  Figures  6,  7,  and  8. 
The  normal  force,  Magnus  force,  and  axial  force  (without  base  drag)  are  virtu¬ 
ally  unchanged  using  constant  stepsizes  of  .008  and  .027  cal.  A  third  case 
shown  in  these  same  figures  used  a  stepsize  of  .027  cal.  at  the  starting  plane 
and  increased  by  5%  every  10th  step  to  .053  cal.  at  the  base  of  the  projec¬ 
tile.  The  previous  PNS  solutions  for  this  case  used  stepsizes  in  the  range 
from  .008  to  .01  cal.  A  considerable  amount  of  computational  time  is  saved 
due  to  the  fewer  number  of  steps  required  to  inarch  completely  over  the  body. 

The  effect  of  smoothing  on  these  same  forces  at  Mach  3  is  demonstrated  in 
Figures  9,  10,  and  11.  Table  1  lists  the  values  of  smoothing  used  to  march 
each  solution.  Definition  of  the  smoothing  parameters  discussed  here  is  made 
in  Reference  12.  Run  #1  represents  the  (minimal)  type  of  smoothing  needed  to 
solve  for  a  body-alone  configuration  at  low  angle  of  attack.  Run  #6  repre¬ 
sents  the  typically  large  amounts  of  smoothing  needed  for  high  angle  of  attack 
solutions  or  finned  body  configurations.  It  is  noteworthy  that  the  normal 
force  is  most  affected  by  the  excessive  amounts  of  dissipation  being  forced 
into  the  solution.  On  the  other  hand,  axial  force  and  Magnus  coefficients 
appear  to  be  less  affected  by  excessive  smoothing.  For  Run  #6,  the  apparent 
lack  of  sensitivity  of  axial  force  to  the  excessive  amounts  of  smoothing  is 
caused  by  an  offsetting  of  the  pressure  and  viscous  components  of  drag.  The 
pressure  drag  experienced  an  increase  of  10%  over  Run  #1,  while  the  viscous 
drag  decreased  by  30%. 


TABLE  1.  MARCHING 

SOLUTIONS; 

COMPARISON 

OF  SMOOTHING  PARAMETERS  S0CBT, 

M 

=  3,  a  = 

2°,  P0/V  = 

.19  (20K  RPM) 

Solution  # 

OX  (cal) 

SMU 

SMU I M  EPSA 

EPSB 

1 

.027 

.005 

0.0  0.0 

0.0 

2 

.027 

.010 

0.0  0.0 

0.0 

3 

.027 

.100 

.25  0.0 

0.0 

4 

.027 

.200 

.50  0.0 

0.0 

5 

.027 

.200 

.50  0.1 

0.2 

6 

.027 

1.000 

2.00  0.1 

0.5 

Figures  12  and  13  compare  the  FINPNS  (Solution  #2)  and  PNS  solutions  to 
experimental  measurements  at  Mach  3.  These  graphs  show  the  development  of  the 
normal  and  Magnus  force  coefficients  over  the  body  and  provide  additional 


12.  D.  S.  Chaussee,  J.  L.  Patterson ,  P.  Kutler ,  T.  Pulliam ,  and  .7.  L.  Steger , 
"A  Numerical  Simulation  of  Hypersonic  Viscous  Pious  Over  Arbitrary 
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validation  of  the  results  from  FINPNS.  In  addition.  Table  2  lists  values  of 

the  components  of  Magnus  force  coefficient;  i.e.,  C_  ,  wall  pressure  comoo- 

Pw 

nent;  C  ,  circumferential  wall  shear;  and  C  ,  longitudinal  wall  shear.  For 

x 

all  cases,  an  adaptive  grid  technique  such  as  that  mentioned  in  Reference  2 
was  used  to  control  the  grid  resolution  at  the  wall. 


TABLE  2. 

COMPARISON  OF 
M  =  3,  a  = 

MAGNUS  FORCE  COMPONENTS  SOCBT, 

2°,  PD/V  =  .19  (20K  RPM) 

it 

El 

FINPNS  (#2) 

+.097  x  10“5 

.213  x  10'3 

-.285  x  10-2 

-.00264 

PNS 

-.147  x  10-5 

.180  x  10‘3 

-.278  x  10"2 

-.00259 

EXPERIMENT 

(Ref.  10) 

-.00270 

Figure  14  compares  axial  force  coefficients  (without  base  drag)  for  PNS 
and  FINPNS  calculations.  The  discrepancy  between  the  two  solutions  was  found 
to  be  caused  by  the  fact  that  the  FINPNS  force  package  makes  use  of  a  second 
order  Taylor  series  expansion  of  the  velocities  at  the  wall  to  define  velocity 
gradients  for  calculations  of  shear  stresses,  while  the  previous  PNS  calcula¬ 
tions  employed  a  first  order  approximation.  The  FINPNS  code  was  re-run  using 
a  first  order  Taylor  series  expansion  to  define  wall  velocity  gradients  and 
results  for  the  viscous  component  of  axial  force  compared  with  the  second 
order  accurate  FINPNS  result  and  with  the  previous  PNS  result,  as  shown  in 
Figure  15.  Both  of  the  first  order  predictions  of  axial  viscous  force  compare 
well,  and  the  difference  between  the  first  and  second  order  prediction 
accounts  for  the  discrepancy  in  the  total  axial  force.  Normal  and  Magnus 
forces,  which  are  primarily  pressure  forces,  were  not  strongly  dependent  on 
the  order  of  the  Taylor  series  expansion  used  for  wall  velocity  gradients. 

Figures  16,  17,  and  18  show  the  pitching  moment  coefficient,  the  pitch- 
plane  center  of  pressure,  and  Magnus  moment  coefficient  as  a  function  of  Mach 
number  as  predicted  by  the  FINPNS  and  PNS  codes,  compared  with  the  experiment. 
Excellent  agreement  between  the  predictions  of  the  two  codes  and  experiment 
are  seen  for  the  pitch-plane  center  of  pressure  and  the  pitching  moment  coef¬ 
ficient.  Good  agreement  is  seen  between  the  predicted  and  experimental  values 
for  the  Magnus  moment  coefficient  at  the  higher  values  of  Mach  number,  though 
this  agreement  falls  off  somewhat  as  the  Mach  number  decreases.  These  results 
further  document  the  suitability  of  the  technique  for  the  prediction  of  the 
supersonic  flow  about  spinning  projectiles. 

B.  Finned  Projectile  Configuration 


The  finned  body  configuration  for  which  calculations  have  been  performed 
resembles  closely  the  M735  Army  projectile.  The  modeled  finned  body  configur¬ 
ation  is  characterized  by  a  conical  nose  section  joined  to  a  smooth 
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cylindrical  main  body  with  six  symmetrical  swept  fins  attached  to  the  aft 
section  of  the  projectile.  Figure  19  displays  the  basic  dimensions  of  this 
configuration.  The  actual  projectile  differs  from  the  modeled  projectile  in 
that  the  actual  projectile  has:  (1)  circumferential  grooves  over  much  of  the 
cylindrical  portion  of  the  body  to  prevent  the  sabot  from  sliding  off  the  body 
in  the  gun  tube;  (2)  fins  which  have  a  non- symmetrical  sectional  geometry  to 
induce  roll;  and  (3)  a  slightly  rounded  nose.  Modeling  the  projectile  with  a 
sharp  nose  and  a  symmetrical  fin  section  are  not  the  result  of  inherent  limi¬ 
tations  of  the  computational  model,  but  rather  a  matter  of  convenience  for 
these  initial  calculations.  While  modeling  of  the  sabot  grooves  may  be  possi¬ 
ble  using  surface  blowing,  wind  tunnel  results  have  shown  that  such  grooves 
have  almost  no  effect  on  the  value  of  normal  force  and  pitching  moment.  These 
grooves  do,  however,  have  a  noticeable  effect  on  drag,  particularly  at  higher 
angles  of  attack. 13 

A  shadowgraph  of  the  actual  projectile  in  flight  at  Mach  4.3  is  shown  in 
Figure  20  and  displays  some  of  the  relevant  features  of  the  flow  field;  a  bow 
shock  wave  emanating  from  the  nose  of  the  projectile,  shocks  at  the  leading 
edge  of  the  fins,  expansion  waves  at  the  cone-cylinder  junction,  and  a  bound¬ 
ary  layer  which  increases  along  the  body. 

Results  are  presented  here  for  Mach  numbers  of  3,  4,  and  5,  two  degrees 
angle  of  attack,  and  turbulent  flow  conditions  over  the  body.  Atmospheric 
flight  conditions  were  simulated  by  maintaining  the  body  temperature  at  the 
free  stream  value  of  294  K.  Calculations  were  made  with  two  of  the  fins 
oriented  vertically,  enabling  a  half  plane  of  symmetry  to  be  applied. 

Generation  of  an  initial  plane  of  data  using  the  conical  starting  proce¬ 
dure  was  performed  at  a  position  .36  calibers  from  the  tip  of  the  cone.  Solu¬ 
tions  were  first  converged  for  large  values  of  explicit  and  implicit  smoothing 
(see  Table  3)  and  typically  required  180  steps.  Implicit  smoothing  was  then 
removed,  explicit  smoothing  reduced  to  a  value  of  0.015,  and  the  solution 
reconverged  with  approximately  an  additional  450  steps.  The  starting  solution 
for  each  Mach  number  was  obtained  using  a  step  size  of  .0072  calibers,  19 
circumferential  points  spaced  equally  in  the  half  plane,  and  45  constantly 
stretched  points  from  the  body  surface  to  the  shock.  The  spacing  from  the 
wall  to  the  first  point  above  the  wall  was  adjusted  so  that  the  first  point 
above  the  wall  was  in  the  laminar  sublayer. 

Using  this  initial  plane  of  data,  the  solution  was  then  marched  down  the 
cone  and  onto  the  cylinder  to  a  position  approximately  one  caliber  in  front  of 
the  fins.  At  this  point  the  circumferential  gridding  was  increased  from  19  to 
121  equally  spaced  points  to  improve  the  resolution  for  marching  over  the 
fins,  and  the  solution  marched  to  the  beginning  of  the  fins.  Over  the  cone- 
cylinder  portion  of  the  body  a  step  size  of  .0072  calibers,  explicit  smoothing 
of  .015,  and  45  points  constantly  stretched  points  from  the  body  to  the  shock 
were  used.  Spacing  from  the  wall  to  the  first  point  above  the  wall  was  again 
maintained  so  that  one  point  was  in  the  laminar  sublayer. 


13.  F.  Brandon ,  "private  communications ,"  U.S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground ,  Maryland. 


TABLE  3.  SMOOTHING  PARAMETERS  FOR  FINNED  BODY  CALCULATIONS 


SMOOTHING 

CONE 

STARTS 

MARCHING  MODE 

PARAMETERS 

1st  Solution 

Final  Solution 

Cone-Cylinder 

Fins 

M  =  3 

M  =  4 

M  =  5 

SMU 

1.0 

0.015 

0.015 

1.0 

1.0 

1.0 

SMUIM 

2.0 

0. 

0. 

2.0 

2.0 

2.0 

EPSA 

1.0 

0. 

0. 

0.1 

1 .-1.5 

1.5 

EPSB 

1.0 

0. 

0. 

0.5 

1.-1.5 

1.5 

Once  on  the  finned  portion  of  the  body,  the  grid  was  obtained  through  the 
use  of  an  elliptic  grid  generator.5  Points  on  the  body  surface  were  clustered 
near  the  leading  edge  of  the  fins,  as  shown  in  Figure  21.  Figures  22a  and  b 
show  a  cross  section  of  the  grid  on  the  finned  portion  of  the  body  at  an  axial 
location  of  X/D  =  13.2.  Grid  points  are  clustered  near  the  body  to  resolve 
the  boundary  layer. 

In  marching  the  solution  over  the  finned  portion  of  the  body,  substantial 
difficulty  was  encountered  in  obtaining  a  solution,  particularly  at  the  axial 
locations  near  the  beginning  of  the  fins  and  at  the  axial  location  where  the 
fins  reach  their  maximum  span.  In  addition  to  adding  significant  amounts  of 
explicit  and  implicit  smoothing  (see  Table  3),  the  distance  from  the  wall  to 
the  first  point  above  the  wall  had  to  be  increased  to  a  value  ten  times  that 
used  immediately  in  front  of  the  fins  in  order  to  obtain  a  solution  for  a  step 
size  of  .0072  calibers.  Decreasing  the  step  size  by  a  factor  of  ten  allowed 
the  radial  spacing  at  the  wall  to  be  reduced  by  a  factor  of  three,  but  the 
resulting  rise  in  computational  time  prohibited  this  approach.  Solutions  were 
obtained  at  Mach  4  and  5  for  lower  levels  of  smoothing,  but  oscillations  of 
the  pressure  on  the  leading  edge  of  the  fins  were  evident.  The  higher  levels 
of  smoothing  used  damped  these  oscillations  out  without  affecting  the  predic¬ 
tion  of  forces  and  moments  significantly. 

Figure  23,  which  displays  the  pressure  along  the  body  on  the  wind  and  lee 
sides  and  at  a  point  in  between,  sheds  some  light  on  the  reasons  for  the  comp¬ 
utational  difficulties  over  the  finned  portion  of  the  body.  A  sharp  rise  in 
pressure  on  the  leading  edge  of  the  wind  and  lee  fins  is  seen,  followed  by  a 
sharp  drop  at  the  axial  location  where  the  fin  reaches  its  maximum  span. 

Figures  24  and  25  show  the  development  of  the  normal  force  coefficient  and 
pitching  moment  coefficient  (referenced  to  the  center  of  gravity  position, 
shown  in  Figure  19)  over  the  body.  The  normal  force  coefficient  shows  a 
moderate  contribution  due  to  the  conical  nose  and  cylindrical  portion  of  the 
body  and  larger  contribution  due  to  the  finned  portion  of  the  body.  The 
development  of  the  pitching  moment  coefficient  over  the  body  demonstrates  the 
stablizing  influence  of  the  fins,  changing  the  pitching  moment  from  a  positive 
value  to  a  large  negative  (stable)  value  at  the  aft  end  of  the  fins.  In  these 
figures  comparison  is  made  with  results  obtained  with  the  NSWC  inviscid  code 


(SWINT}1**  and  LCWSL  data15  base  values  for  the  actual  projectile  configura¬ 
tion. 

The  SWINT  code  -  SWINT  is  an  acronym  for  Supersonic  Wing  INlet  Tail  - 
solves  explicitly  the  Euler  equations  for  supersonic  flow  over  bodies  with 
fins  and/or  wings,  and  for  the  external  flow  about  bodies  with  inlets.  The 

code  makes  use  of  the  thin  fin  approximation,  collapsing  each  fin  along  a 

single  radial  plane  in  the  grid.  Fin  thickness  can  be  accounted  for  in  the 
code  by  application  of  the  appropriate  local  analysis  such  as  shock  compres¬ 
sion  and  Prandtl-Meyer  expansion  theories.  The  fin  edges  must  be  sharp  and 

cannot  extend  beyond  the  bow  shock.  Additionally,  the  flowfield  must  remain 
supersonic  throughout  the  entire  computation. 

Development  of  the  normal  force  and  pitching  moment  coefficients  as  pre¬ 
dicted  by  the  FINPNS  and  SWINT  codes  compare  well.  Good  agreement  between  the 
total  value  of  normal  force  and  pitching  moment  for  both  procedures  is  seen 
compared  with  the  LCWSL  data. 

Pressure  distributions  on  the  wind  and  lee  sides  of  the  $  =  120°  fin  at 
the  1/4  and  1/2  span  positions  are  shown  in  Figure  26 a  and  b  compared  with  the 
SWINT  code  predictions.  Differences  in  the  predicted  pressures  at  the  leading 
edge  of  the  fin  are  attributable  to  the  thin-fin  approximation  applied  for  the 
SWINT  calculations.  The  FINPNS  code  seems  to  predict  a  more  even  distribution 
of  lift  along  the  chord  of  the  fin,  while  the  SWINT  predictions  show  slightly 
more  lift  being  developed  at  the  leading  edge  of  the  fin. 

The  variation  of  the  predicted  values  of  the  slope  of  the  normal  force  and 
pitching  moment  coefficients  with  Mach  number  is  compared  with  the  values  pre¬ 
dicted  by  the  SWINT  code  and  with  LCWSL  data  in  Figures  27-28.  The  FINPNS 
code  is  seen  to  overpredict  the  normal  force  and  pitching  moment  coefficients 
compared  with  the  LCWSL  data,  while  the  SWINT  code  gives  generally  better 
agreement. 

The  variation  in  the  predicted  center  of  pressure  with  Mach  number  is 
shown  in  Figure  29  and  is  also  compared  with  values  predicted  by  the  SWINT 
code  and  with  LCWSL  data.  Both  codes  predict  similiar  values  of  center  of 
pressure  and  show  good  agreement  with  the  LCWSL  data. 


IV.  CONCLUSIONS 

The  application  of  the  Parabolized  Navier-Stokes  Technique  to  a  Secant- 
Ogi ve-Cy 1 i nder-Boattai 1  (SOCBT)  projectile  has  verified  recent  code  improve¬ 
ments  which  allow  the  code  to  be  utilized  in  an  even  more  computationally 
efficient  manner.  This  report  has  demonstrated  that  consistently  accurate 


14.  A.  R.  Wardlaw ,  Jr.,  F.  P.  Baltakis,  J.  m  Soloman,  and  L.  B .  Hackerman , 

"An  Tnviscid  Computational  Method  for  Tactical  Missile  Configurations 
NSWC  TR  81-457. 

15.  Unpublished  range  lata,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland. 
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results  can  be  achieved  for  the  SOCBT  configuration  at  small  angles  of  attack 
over  a  range  of  supersonic  Mach  numbers,  establishing  a  satisfactory  compari¬ 
son  against  benchmark  results. 

Results  of  the  FINPNS  calculations  for  the  finned  projectile  configuration 
have  shown  fair  agreement  with  LCWSL  data  and  Euler  (SWINT)  code  predictions 
but  further  improvements  are  required  before  a  satisfactory  predictive  capabi¬ 
lity  can  be  considered  to  exist.  Current  levels  of  smoothing  should  probably 
be  reduced  and  better  resolution  of  the  boundary  layer  on  the  fins  is 
required.  Implementation  of  a  global  iteration  technique  may  help  overcome 
some  of  these  difficulties. 
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initial  data  plane 


AU  DIMENSIONS  M  CALBERS  (ONE  CALIBER  =  57.2  mm) 

Figure  2.  Secant  Ogive-Cylinder-Boattail  (SOCBT)  Configuration 
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Circumferential  Surface  Pressure  Distribution  on  SOCBT 
Configuration;  M  =  3,  a  =  2°,  PD/V  =  0,  Rep  =  1.2  million 
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Figure  7.  Total  Magnus  Force  Coefficient  Developing  over  SOCBT  Configuration; 

M  =  3,  a  =  2°,  PD /V  =  .19,  ReD  -  1.2  million 
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Figure  8.  Total  Axial  Force  Coefficient  Developing  over  SOCBT  Configuration; 

M  =  3,  a  =  2°,  PD/V  =  .19,  Ren  =  1.2  million 
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Figure  9.  Total  Normal  Force  Coefficient  Developing  over  SOCBT  Configuration 

M  =  3,  a  =  2°,  PD /V  =  .19,  ReD  =  1.2  million 


Figure  10.  Total  Magnus  Force  Coefficient  Developing  over  SOCBT 

Configuration;  M  =  3,  a  =  2°,  PD/V  =  .19,  Rep  =  1.2  million 
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Figure  11. 


Total  Axial  Force  Coefficient  Developing  over  SOCBT  Configuration; 
M  =  3,  a  =  2°,  PD/V  =  .19,  ReD  =  1.2  million 
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Figure  12.  Total  Normal  Force  Coefficient  Developing  over  SOCBT 

Configuration;  M  =  3,  a  =  2°,  PD/V  =  .19,  Rep  =  1.2  million 
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Figure  13.  Total  Magnus  Force  Coefficient  Developing  over  SOCBT 

Configuration;  M  =  3,  a  =  2°,  PD/V  =  .19,  Rep  =  1.2  million 


Figure  14.  Total  Axial  Force  Coefficient  Developing  over  SOCBT  Configuration; 

M  =  3,  a  =  2°,  PD/V  =  .19,  Rep  =  1.2  million 
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Figure  15.  Viscous  Component  of  Drag  Coefficient  Developing  over  SOCBT 
Configuration;  M  =  3,  a  =  2°,  PD/V  =  .19,  Reg  =  1.2  million 


Figure  16.  Slope  of  Pitching  Moment  Coefficient  versus 
Mach  Number  for  SOCBT  Configuration 


Figure  23.  Axial  Pressure  Distribution  at  4>  =  0°,  4>  =  90°,  and  4>  = 

M  =  4,  a  =  2°,  Rep  =  3.2  million 
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Figure  24.  Development  of  Normal  Force  Coefficient  along  Body 

M  =  4,  a  =  2°,  Ren  =  3.2  million 
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Figure  26b.  Chordwlse  Pressure  Distribution  at  1/2  Span  Position 
4  -  120°  Fin,  M  =  4,  a  =  2°,  ReD  =  3.2  million 
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Figure  27.  Zero-Degree  Slope  of  the  Normal  Force  Coefficient 

versus  Mach  Number 


Figure  28.  Zero-Degree  Slope  of  the  Pitching  Moment  Coefficient 

versus  Mach  Number 
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Figure  29.  Normal  Force  Center  of  Pressure  versus  Mach  Number 
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LIST  OF  SYMBOLS 


speed  of  sound 

axial  force  coefficient 

wall  shear  component  of  axial  force  coefficient 

normal  force  center  of  pressure 

pitching  moment  coefficient 

dCm/da,  slope  of  the  pitching  moment  coefficent 

Cn/(PD/V)a,  slope  of  the  Magnus  moment  coefficient 

Magnus  (yawing)  moment  coefficient 
normal  force  coefficient 

dC^/da,  slope  of  the  normal  force  coefficient 
surface  pressure  component  of  Magnus  force  coefficient 
Magnus  (side)  force  coefficient 

longitudinal  wall  shear  component  of  Magnus  force  coefficient 
circumferential  wall  shear  component  of  Magnus  force  coefficient 
diameter  of  model 

total  energy  per  unit  volume  of  fluid,  normalized  by  p^a* 
flux  vectors  of  transformed  gasdynamic  equation 
projectile  length 
Mach  number 

pressure  normalized  by  p.a* 
non-dimensional  spin  rate  about  model  axis 
Reynold's  number  based  on  diameter,  pJM^aJD/Ue, 
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A 

Re 

A 

s 


u,v,w 


x.y.z 

a 

U 

S.n.c 


p 

♦ 


Reynolds  number,  p^ajl/w^ 
viscous  flux  vector 

Cartesian  velocity  components  along  the  x,  y,  z  axis, 
respectively,  normalized  by  a„ 

physical  Cartesian  coordinates 

angle  of  attack 

coefficient  of  viscosity,  normalized  by  free  stream  value, 

computational  coordinates  In  the  axial,  circumferential,  and 
radial  directions 

density,  normalized  by  free-stream  density 
circumferential  angular  coordinate 


Subscripts 

•  =  free-stream  conditions 

wall  *  body  surface  values 
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